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We define different classes of local random quantum circuits (L-RQC) and show that: a) statistical 
properties of L-RQC are encoded into an associated family of completely positive maps and b) 
average purity dynamics can be described by the action of these maps on operator algebras of 
permutations (swap algebras). An exactly solvable one-dimensional case is analyzed to illustrate 
the power of the swap algebra formalism. More in general, we prove short time area-law bounds 
on average purity for uncorrelated L-RQC and infinite time results for both the uncorrelated and 
£f) ■ correlated cases. 
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The quantum-mechanical state space of a many-body system (say N particles) is huge, as its dimension scales 
exponentially with N. The exact specification of any quantum state therefore requires, its face value, an exponential 
■ number of parameters, a daunting task both numerically and experimentally. As in classical statistical mechanics one 
is then naturally led to resort to probabilistic methods to explore the typical properties of ensembles of quantum states. 
The question now becomes: what are the "natural" probability distribution laws one should consider over the state 
space? A first conceptually and mathematically appealing answer is: the group-theoretic Haar distribution This 
latter choice in fact seems to implement the principle of minimal a priori assumptions and at the same time allows one 
to exploit the powerful computational tools of representation theory [l[ . On the other hand the Haar measure presents 
problems from a quantum information theoretic perspective. In fact it has been known for almost a decade now that 
^ ■ sampling exactly the Haar distribution over a many-body Hilbert space (say N particles) by a quantum circuit made 
out of local gates i.e., a quantum computer, is an exponentially (in N) hard task Q. This result has been recently 
strengthened by considering all possible states that are achievable by time-dependent local Hamiltonians acting for a 
^C) ' polynomial amount of time: the quantum states that are produced in this way occupy an exponentially small volume 
, in Hilbert space Q. This implies that the overwhelming majority of states in Hilbert space are not physical as they 
' can only be produced after an exponentially long time. This startling fact led the authors of Q to dub the Hilbert 
space " a convenient illusion" (sic). 

However, generic properties of Haar distributed quantum states play a crucial role in the recent revival of interest 
on rigorous foundations of statistical mechanics started with the seminal papers [H and @ . In the light of the above 
mentioned findings about the unphysical nature of Haar distributed random states, this foundation for statistical 
ensembles appears to be conceptually questionable. How does one introduce a family of "physically accessible" states 
endowed with a natural probability measure? In Ref. @ a first attempt to undertake this endeavor was made by 
. . ' introducing a family of random matrix product states, in this paper we will focus primarily on local random quantum 
circuits (L-QRC) families 0413 ■ Roughly speaking (a precise mathematical definition will be given in the next 
section) a L-RQC is a random sequence of k random unitary transformations (quantum gates) each of which acts on 
small number of localized quantum degrees of freedom. The states of N particles obtained by the action of a L-RQC 
(with k = Poly(iV) ) on an initial e.g., factorized, fiducial one are now regarded as the physically accessible ones. 

A second important physical motivation for the introduction of L-RQC comes from the perspective of quantum 
control and quantum information processing with limited resources. In this case lack of ideal experimental ability 
may prevent one to know the precise structure of the quantum gate one has enacted as well, in the case of bounded 
spatial resolution, the precise set of local degrees of freedom e.g., qubits involved. This scenario naturally leads to a 
modelization in which both unitary transformations and their supports are effectively random variables. 

A third compelling physical motivation behind the introduction of L-RQCs is that they provide a natural way to 
"simulate" the quantum dynamics of local random (time-dependent) Hamiltonians with a discrete circuit. This, for 
example, can be seen by expanding the evolution generated by a local-Hamiltonian using a Suzuki- Trotter type of 
factorization Q . In this way the study of the ensembles for circuits whose length increase unboundedly i.e., the limit 
k — > oo, will shed light on the fundamental problem of quantum equilibration [14j. Finally, L-RQC have been recently 
considered in the context of t— designs theory @, §]■ Even if this latter one is a quite important area of application of 
L-RQCs, t— design theory is not the focus of the present work; it will briefly touch upon just in the last section. 

The goal of this paper is to provide a mathematically precise account of the physical approach to L-RQC outlined in 
Refs. [3. [Tot- Generalizations of the results contained therein and rigorous proofs of a few conjectures are presented. 
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II. LOCAL RANDOM QUANTUM CIRCUITS 

Let us start by laying down the mathematical framework we are going to elaborate on. Let {/ii}ieA be a collection 
of "local" Hilbert spaces labeled by a set A with finite cardinality |A|. For each subset C A we have an affiliated 
Hilbert space T-Lq := <E>i£nhi, as well as the operator algebra A(£l) := B(Hq). If fli C fl 2 we will consider A(£li) C 
A(ft 2 ) by identifying a x e A{£li) with a 2 = ai ® ln 2 -ni G A(ft 2 ). In particular A(Q) C A(A), (Vf2 C A) and 
.4(A) = [J ncA A(fl). For simplicity we will further assume that each local space has the same finite dimension d i.e., 

hi = C d , (Vi € A). In this case the global state-space Ha is isomorphic to (C d )® ' A ' (dim Ha = d' A '). 
Definition 1 

A local random quantum circuit family Cfc[S, C, q^ k '] of length k is defined by the following data: 

i) A family of probability densities S := {d/io}ncA over the groups U(il) := {U S A(£l) / W = U^ 1 } 

ii) A subset C of the power set 2 A . Elements of £ will be referred to as the local regions. 

ii) A probability law qW : C k -> [0, 1] / £1 := (0 Xj . . . , H fc ) ->■ ? (fc) (^), Ene£« = 1- 

A local quantum circuit of length A: is a unitary in U(A) of the form Uci '■= UkUk-i ■ ■ - U\ where Vz G {1, ...,&}=> 
a)[/i € M(fi<); b) O € £ fe . 

Roughly speaking, the idea is that the C/o's are unitary- valued random variables distributed according the law 
q( k \fl)dUn (dllfi = JliLi ^Mn^^sii))- Two different layers of randomness are involved in our construction of L-RQC, 
they are associated with items i) and ii) in Def. 1 respectively. The first one is a stochastic process (of length k) in 
which k local regions of the base set A are selected according to the joint probability distribution q( k K The second 
is the selection, according to the d/zo's of k unitaries acting on the state-space of each of the local regions randomly 
selected in the previous step. 

The physical picture behind these definitions is quite simple: the experimenter has access only to random unitaries 
localized over regions ft belonging to a distinguished subset £ C 2 A . The definition of the local regions set C is critical 
and it is where extra (physically motivated) input is needed. The idea is that the f2 £ C are "small local " subsets 
of A. In the following we will define RQCs over (hyper) graphs, in that case A is the vertex set and £ will coincide 
with the set of (hyper)edges; in general one may think of the local regions as sets of cardinality O(l) (where |A| is 
the large scaling quantity). Also, in the case in which A is naturally equipped with a metric structure e.g., the graph 
theoretic distance, the local regions will also be required to have small diameter i.e., diam(fi) = 0(1). 

In the following we will almost exclusively restrict ourselves to the case where all the d/xn's are the Haar measure 
over U(Vl). This symmetry assumption is crucial for our analysis as it allows one to resort systematically to powerful 
group representation theory tools [l|. Accordingly we will hereafter drop the 5 symbol in the definition of Cfe. Notice 
however that other choices are possible. For example in [l3| the d^s were concentrated on a finite set of universal 
gates. A first goal for this general approach is to identify families of L-RQC that are at the same time physically 
motivated and amenable to rigorous mathematical analysis. A couple of L-RQC families that stands out as a natural 
target for further investigations is given by the following types of 

Definition 2 

• Time dependent Markovian L- RQC: 

Q«(n) = M( fe - 1 )(o fe ,n fe _ 1 )M( fe - 2 )(o fc _ 1 ,n fe _ 2 )---M«(Q 2 ,n 1 ) g «(n 1 ) 

where {M^} k r\ are k — 1 stochastic |£| x \£\ matrices and q( ' is a distribution over C. If = Af, (Vj) we 
have a time- independent Markov process over C 

• Uncorrelated L-RQC: q^(ft) = Ilj=iP^H^j)i where the are k distribution laws over C. If q^ — q, Vj we 
have a time-independent uncorrelated process. 

Even in these restricted setups the landscape of possibilities is huge and one has to resort, once again, to physical 
input to further specialize the models. For example if A is the vertex set of a graph Q := (A, E) and C = E one may 
consider such that Mw (f2i, f2 2 ) ^ iff the edges Oj and fl 2 share at least one vertex. In this way the AfW's 

define a sort of random walk on the graph Q. 
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In passing we notice that one can use Ck[C, q( k >] to turn A(A) in to a noncommutative probability space with 
the expectation <f> : A{A) — > C / A — > (j)(A), where 




Here (A, B) :— Tr (A^B^j the standard Hilbert-Schmidt scalar product over *4(A) and uj is an distinguished (initial) 
density operator in „4(A). 

A physically compelling model of L-RQC that received much attention recently is based of graph structures Q = 
(A, E) |7H9l. Il3|. In our formalism this amount to say that i) A is the set of vertices of a graph Q, ii) C = E C A 2 is 
the set of edges of Q. Most of the literature on random quantum circuits so far has focused on this graph case in the 
time-independent uncorrclatcd case 2) and with Q being the complete graph. One-dimensional correlated cases have 
been discussed in [nj • 



III. ENSEMBLE MAPS 



In this section we will introduce and start to analyze a set of completely positive (CP) maps 16] naturally affiliated 
with any L-RQC family 0, [HI- These maps comprise the full statistical content of the L-RQC family and play a 
central role in the following of the paper. 

Using a compact notation we can write ([1]) in the symbolic form 4>(A) = (u, Ut^U) . For any fixed observable 
A G A(A) the mapping Xa ■ XJ (to, XJ^AXJ) defines a classical (commutative) random variable. The statistics 
of Xa describes the fluctuations of the (purely) quantum expectation (a;, LJtylU) over the ensemble of the local 
random circuits U in Ck[C,q^]. While the p-moments of the noncommutative random variables A (p € N) are 

defined by n P (A) — <fi(A p ) the moments of the corresponding classical variables are given by h p (Xa) = Xa{V) p = 

(lj, UtvlU) p = (w® p , (TJ1'AU)®p ). These moments can be conveniently expressed in terms of a family of maps 1Z P 

associated with C k [C,q^} : n P (X A ) := (uj® p , K p {A® p )) U^M- Here K p : A(A)®p -> A{A)® P / A -> (fUp" 
are completely positive (CP) maps [la ]. More explicitly, if A g A(A)® P , one can write 

Tl p (A)= ]T q^(n)TZ p ,n(A), n p ,n(A):= [ dUn(U^f p AU^ (2) 
nec k 

We will refer to the lZ p s as the ensemble maps. Indeed, the distribution law for any A^'s is determined by 
the Fourier transform of the characteristic function XA{t) '■— ^pT^pi-^-A.)- This can be in turn expressed as 

(woo, ftoo(*)(Ax>)) where x x := ©£L :r® p , (x = A, u) and TZoo(t) := ®^L ^TZ P is a formal CP-map over the full 
(operator) Fock space Aoo(A) = ©£L .A(A)® 3 '. 

The ensemble maps TZ p s are clearly uniquely defined by the data in Ck[C,q^] and should be characterized for 
physically relevant L-RQC families. In this paper we will often consider the 1Z P as operators over the Hilbert-Schmidt 
space A(n)® p ; their norms are defined accordingly i.e., \\K P \\ := swp{\\K p (A)\\ / \\A\\ := VTrA^A = 1, A € A{Q.)}. 
Here below we give a summary of their general properties. 



Proposition 

i) The maps 7Z P are: trace-preserving (tr7vL p (X) = tr(X)) and unital (lZ p (t) = 1). 

ii) If g( fc '(f2) = gw(O), (ft := (fife, . . . , fli)), Vfi, the TZ p s are hermitean with respect to the Hilbert-Schmidt scalar 
product over .4(A)® p . 

iii) In the uncorrelated case (Def. 2) they factorize: K p = Hj =1 *R% ) where K p j) (A) = 

Enec1 U) (m P ,n(A), U p .n(A) := J dU n (U&® p AU® p . 

iv) If dUu is the Haar measure and that the distinguished state w is a pure state completely factorized over A i.e., 
w = ®ieA|0i)v&|. Then 

(^ p ,n p , n (A)) = {PW,A) (3) 

where ojq = ®nzn\4>i){4 > i\ an d P p + ^ — So-es a ' s ^ ne (normalized) projector onto the totally <S p -symmetric 

subspace of H® p . 
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v) In the Haar measure case the TV P are positive semi-definite operator ({X, 1Zp (X)) > 0,VA), their spectra are 
contained [0, 1] and always contain 1. 

Proof.— i)~iii) follow easily from Eq ©. 

iv) In this case representation theory implies that the 7?. Pj n's (in [21 k — 1) are projections on the commutant of 
the representations U G U(Jln) — > U® p G U(Hq P ) [l|. By the Schur-Weyl duality this commutant is the algebra 
generated by the permutations S p 3 a : — > 'Hq P / ®f = i <i>i ~ * ®i=i0o-(i) 0- Moreover, in the computations 
of the moments one has eventually to contract with ui® p whose support is entirely contained in the iS p -symmetric 
subspace of This implies that the only relevant part of the projection is the one associated with the identity 
irrep of S p . This leads to the explicit formula ([3]). 

v) In the Haar measure case each of the IZp^ in hi) of Prop. is a convex combination of the projectors lZ Pt n 
and it is therefore a positive semi-definite. Moreover , \\TZp , || = \\ YlneC 5^'K^)^p>oll — J2nec I^K^W^-p^W — 
T,nec = h hence Sp(^p } ) C [0, 1]. Since e.g, 1 is always a fixed point 1 G Sp(^ j) ). □ 



IV. FIRST MOMENTS 



In this section we will analyze the L-RQC first moments of quantum observables by means of the ensemble maps 
1Z±. In particular the dynamics of these first moments can be studied in the limit in which the circuit size grows 
unboundedly i.e., k oo. 

Let to G -4(A) be a density matrix (i.e., u > 0, Tr(w) = 1) and A G -4(A) a quantum observable. The expectation 
value of the random variable U i— > (w, AU) over the ensemble of uncorr elated, time-independent circuits of length 
k is given by (see previous Sect. ) by 4>k{A) := (w, 7£f (A)), where IZi — J2nec 1 0^)"^i,n< We will also assume that 
£ is a covering of A i.e., Ung^Jl = A. 



Proposition 1 

i) In the Haar measure case the ensemble maps 7^-1^ are projections: TZi t n(X) = d~' a 'tn ® Trn(X). 

ii) ^i,ni^i,n 2 = "^i,n 2 ^i,Oi =^i,niun 2 - 

iii) Sp(^i) C {q(a) := ^ q^(Q k )a k /ax,..., a [C \ ^2 = {0, 1}} 3 1 

iv) The following estimate for the convergence <j>k{A) — > 4>oo{A) = cM A 'TrA, holds 



k > 



log(||o;|| 2 ||A|| 2 /e) + (|£|-l)log2 



!og- 



\MA) -<f>oo(A)\ < e 



(A) 



where 57, = argmin^ (^)- 
Proof. 

i) Follows from the Schur Lemma [l|. In particular if X = Xq_ ® Xjjc g A(il) (S) -4(f2 c ) one has TZi i n(X) = 
d-l"lln®XnoTr(Xn). 

ii) When fii fl f2 2 = the claim is obvious. Let us assume then fii n fi 2 7^ and w.l.o.g. that A = f?i U fi 2 . We write 
A = Ai U A 2 U A 3 where Ai := fix - fi 2 , A 2 = fli n fi 2 , and A 3 := fi 2 - Oi- Let us consider X G -4(A) = <g>f =1 ,A(Aj) 
of the form X x ® X 2 ® X 3 . One has TZ h nA x ) = dA n ^l Al ® l Aa ® X 3 TrXiTrX 2 , therefore ftx,fi 2 ("£x,Qi PO) = 
dA^\dA^\t Al ^l A2 ^t A3 Til A2 Th:X 3 TkXiTrX 2 ,= dA Q ^-\^\+\^nsi 2 \ 1a jja^ T T x t = dA A \l A T T X, (|A| = |fiiU 
fi 2 | = |fii| + |fi 2 | — \flx n fi 2 |. Computing TZi v n 1 (7£x,fi 2 (X)) clearly gives the same result. 

iii) Since IZi is a convex combination of commuting projections its diagonalization is formally immediate. Let 



us introduce the notations P, 



(0) 



1 



Kn k , and 

(&k 



identity CP map as follows 1 = HtiCo =0 1 ^fc 



= TZn k , ({fii, ■ • ■ >fi|£|} = £)• One can then decompose the 
) = E au ..., a]c ^iU[%H ak) - Therefore 



Ql ,...,a K |=0,l \fc=X / fc=l aeZ l£| 



(5) 



^ j P^'' are necessarily 



where we used 72.i.n fc P[a] = a/c-P[a](Va G ). Notice that not all the projectors P[a] := \\ 
non zero 23], in any case a := (a±, . . . , a Q|C| ) 7^ a' f"[a]-P[a'] = 0. This show that the decomposition is a 
spectral one from which the claim Sp(7?-i) C {<z(a)} aeZ i^i follows. 
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Since 3X such that P[(l, . . . , 1)](X) = TTfeli ^i,n* W = ^ , o (*) = ^I.a(X) = d-l A IH A TrX ^ one sees 

' fc = 1 ^ 

that P[(l, . . . ,1)} 7^ and therefore . . . , 1) = Y^n 1 (^) = 1 always belong to the spectrum of IZi. Moreover 
(Prop. 0, v) 1 the maximum eigenvalue and q(a) = 1 — J^i a-=o 1 0^») < 1 if a 7^ (1, . . . , 1). 

iv) We first observe that the second largest q(a) is upper bounded by 1 — q^(Cl^) where Cl* — argminn q^\Cl). Sec- 
ond, Tli (A) = <?(a) fc P[a](A) = P[(l . . . , 1)}(A) + Ea^(i,...,i) ?(«)^[«](^) &om which lim^ ft* = P[l, . . . , 1] 
and 

- 000(4)1 = |^(4)-(a;,P[l,...,l](4))|< ^ g(a) fc Kw,Pfe](A))| 

Q#(l,...,l) 

< (l- g «(fi*)) fe 2 ||P[a]||||u;||2||4|| 2 <(l-QW(0»)) fe 2l £ l- 1 ||u;|| 2 ||A|| 2 .<e (6) 

Here we used Cauchy-Schwarz |(w,P[a](A))| < ||w|| 2 ||P[a](,4)|| 2 < ||w|| 2 ||P[a]||A|| 2 < ||cj|| 2 ||A|| 2 as P[a] is a projector 
i.e., norm one . Solving the last inequality in (|6|) for k gives the estimate (|4]). □ 

Whereas the asymptotic value 4>oa(A) does not depend on the distribution (as long as U^ec^ = A) the 
convergence rate may well do so. In order have a faster convergence one would like to make #(0*) as large 
as possible i.e., g(0«) = 1/|£|- In this case a weaker form of (UJ) gives k » log(||c<j|| 2 || A|j 2 2l £ l _1 /e) = 

0(|£|log(||w|| 2 ||A|| 2 /e) + |£| 2 log2) -> - </>oo(4)| < e. In several relevant examples one has |£| = 0(|A|) 

so that (jj) gives, for ||4|| 2 = O(l), a convergence estimate polynomial in the system size. 

Of course Q is a very crude estimate that may be strongly improved yielding faster convergence rates. For example 
if A G A(Cl) then P[a]A = for all the a's such that at = when Cl^ C Cl c . This constraint greatly reduces the 
number of non vanishing terms in giving rising to a better bound. In fact, If |{f2 e£/!lc O c }| = O ( | £~2 C | ) then 
|£| in becomes |£| - 0(\Cl c \) = 0(\fl\). 

V. PURITY DYNAMICS AND SWAP ALGEBRAS 

In this section we will show how the ensemble maps 1Z 2 allow one to study the dynamics of average purity and 
how this can be done in terms of dynamically closed subalgebras of permutations. Uncorrelated as well as correlated 
L-RQC will be considered. 

Given a density matrix u> G -4(A) the reduced state wq 6 A(Ct) associated with the region Cl C A is given by 
the partial trace ujq = trn=(w) (Cl c := A — Cl). A family of non-linear functions that plays an important role in 
quantum information theory is given by the a-Renyi entropies: S a (uin) := — — f^- log Tr(w^), (a € R^)- In particular 
lim a _).i S a (oJn) = — Tr (wnlogwn) is the von Neumann entropy that for pure uj quantifies quantum entanglement 
between the two regions CI and Cl c [l7j]. For a = p E N one can write Tr(cj^) = Tr (t^u]® p ^) where : H® p -> 

T-L^ p is the cyclic permutation n : i — > i — 1 modp acting on the Hq factors of H® p = 'Hq P <S>'H® p n - For p — 2 we will 
call this operator the swap associated with CI and it will be denoted by Tq, we will also write 7Z 2j si simply as TZ^i and 
so on. Since = one has Tq = where the Tj act as tt in hf 2 . The case p = 2 has special relevance as 

the 2-Renyi entropy S 2 (l)q) = — logtr(w^) is a convenient way to quantify quantum entanglement that is receiving a 
constantly growing attention from the quantum information and theoretical condensed matter communities see e.g., 
(l8j . If the initial state oj is evolved by a RQC U in Cfc[£, q^] the corresponding purity becomes a random variable 

whose expectation value P k := Tr ((U® 2 )tT U® 2 2 ) U is given by: 

P fe := ^ g«(fi)<u;® 2 , K n (T a )) = (w® 2 , ft(r n )) (7) 
ctec k 

By convexity it follows that S 2 V > — logP^, namely the negative of the logarithm of (J7J provides a lower bound to 
the average entanglement generated by evolving w with RQC in U in Cfc[£, q^ k ']. Two of the key technical insights in 
,9j on which we would like to build upon in this paper are contained in the next two propositions: 

Proposition 2 

The swaps Tq form an abelian group 7a of order 2' A ' whose elements have degree two i.e., Tj 2 = 1. 

T ni Tn 2 = Tq iA q 2 , ttiAft 2 := (Q r - Cl 2 ) U (Cl 2 - Cl x ) {VCl u Cl 2 c A) (8) 
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The group algebra C7a is 2 ' A ' -dimensional abelian subalgebra of C7a C A(A)® 2 under multiplication. C7a will 
referred to as the swap algebra of A. 

Proof.- T ni T n2 = (n i60l T t ) (n je n 2 T,) = (Il^-n, ^) (lL e n in n 2 ??) (ILen,-^ T *) = 
(nie(ni-n 2 )u(n 2 -ni) ^y := ^niAn 2 - Where we have used T 2 — 1, Vi € A. Notice also that the Th's are linearly 
independent in .4(A)® 2 as each of them amounts to a permutation of the product state basis of %® 2 = ®i^{Jif 2 . □ 

Eq. (|SJ) shows that 7a is isomorphic to the power set of A endowed with the internal operation of symmetric set 
difference A. The isomorphism being given by Q — > Tq. Moreover, the swap algebra C7a is isomorphic, as a vector 
space, to the space (C 2 )® ' A ' associated with |A| qubits (Tq = ®jeA|xn(i)), xn characteristic function of O c A ). 
Next proposition shows that C7a is invariant under the 1Zqs. 

Proposition 3 

i) If Q C A let us define its boundary (in £) as dVl := {Q! e C / fT n ft ^ A ft' n tt c ^ 0}, then 

a) Oi £ dn ^> KnATn) = T n , otherwise b) 

UnATn) = a + Xh_ fil + a-PouQi, (9) 

where a± = (c+ ± c~)/2, c± = {d A ± d 5 )/^ 5 ± 1), A = |fii - fi|, B = |fi n 

ii) The swap algebra is invariant under all the ensemble maps 1Z. 

iii) Let P G L(CTa) s.t. T n h> The, (0 c A) then [7e ni , P] = 0, (Vfii C A) 

Proof.— i) a) If Sli C fl c then the result follows trivially from (Ijii ) = Iqi- If fii C f2 then it follows from 
< 2 TnU^ = T f2 _ t2l < 2 T f2l [/ f t 2 = Tn-a.Tn, (r 0l < 2 T ni ) C/^ 2 = T < 2 C/® 2 = T„. 

b) First, one has that (Ih) = Th-ni^ni (JanaJ- In order to compute (TfinnJ, notice that 72.fji (restricted 
on Hq 2 ) is a projection on the 2-dimensional algebra spanned by 1 and Tq x A Hilbert- Schmidt orthonormal basis 
of this space is given by the elements 

F *±T ni 

± ' v /2#^npnrr±T)' 

Hence T^n^Thnni) = E Q =± ^Tr (PaThnoJ ■ Using © one obtains, Tr(F„Tnnni) ~ Tr (Tnnfii + aTm-fi) = 
^lonOil^lOi-oi + Q!d |fi 1 -n| d 2|nnn 1 | j Now the claim f n ows by noticing that |Oi| = |fi a - Q| + |fi H n x | and that, 
from Eq. ©, one has Tq^q 1 Tq 1 = Tnuni- 

ii) It is immediate by observing that IZn = TZfi 1 o 7?.q, o • • • o 7?.Q fc (remind that Jl = (f^, . . . , € £ fe ) and using 
© and ©. 

iii) First notice that (from i)) cr^Q,^) = o^(n c ,Qi), whence Kn^PTn) = ft ni (T n c) = a+(O c , ni)T n . n nj + 

a~(fi c , fiiJTn.^ = a-(n, 1 )T (nufJl) c +a+(0, ni)T (nn n ? }« = a _ (fi, ^i)PPhu^ +a+(^, fii)PTnnn S = Pftm (T n ). 
□ 

Proposition 3 is an important result. It shows that (average) purity dynamics can be mapped onto a dynamical 
problem in a space comprising just |A| qubits. This mapping entails a dimensional reduction (<i 2 )' A ' H >• 2' A L Exploiting 
symmetries this reduction can be made, in some cases, even stronger. To start with, map P in iii) is an involution 
i.e., P 2 = 1, and therefore CTa breaks up in two orthogonal (|A| — l)-qubits subspaces (corresponding to eigenvalues 
±1 of P) that are invariant under all ensemble maps. Moreover, the action of the ensemble maps on the swap algebra 
can be highly reducible thus lowering even further the dimensionality of the purity dynamics problem. References Q 
and [To| contain severale examples (see two below) in which the entanglement dynamics is approached and solved at 
the swap algebra level because of the dramatic decrease of computational complexity. 

A. Uncorrelated case 

In the uncorrelated case the maps TZ^ define a (discrete time) dynamical system on C7a : the element T is mapped 
at "time"= k e N into T(k) := TZ^ oT^ 1 ' • ••ftW(T). In particular in the time independent case T(k) = 1Z k (T). In 
this section we will focus on this latter case. Since the eigenvalues of 1Z are in [0, 1] (see v) in Prop. 0) the asymptotic 
behavior at k = oo is controlled by the eigenspace of 1Z with eigenvalue 1 i.e., Fix (72.) := {X G -4(A) / TZ(X) = X}. 
Next proposition describes the structure of this subspace and gives an explicit formula for the limit purity. 
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Proposition 4 

Let us consider C C 2 A as an hypergraph in A and let {Ci}f =1 be the family of its connected components. [24| i) 



K 

Fix(P) = P| Fix(ftn) = ((g) C{1, T C J) ® A((U? =1 Ci) c ) (10) 
oe£ i=i 

ii) Using the notation of Eq. (J7J, for a totally factorized initial state oj 

K jICI-IO,! , 

Poo- lim(^ 2 ,^ fc (Tn)) = TT ^7H-^ . (^ := O H i = 1, . . . , A) (11) 

i=i 

Proof.- i) Let us first prove the first equality in (p~U|). IfX e Fix(ft) then ||AT|| = ||ft(X)|| = || J2nec <7(^)^P0II < 

En 6 r?(n)ll^nWII < En e£ 9(fi)||ftn||||X]| = (Eo e£ 9(^)) 11*11 = 11*11; Here we used \\K n \\ = 1, (Vfi) and the 
normalization of q. This shows that all these inequalities are actually equalities that in turns implies that, ||72.n(X)|| = 
||Z||,(Vn), and, since the Tin's are projections, Kq{X) = X, (Vfi G £). Therefore Fix(P.) C finer Fix C^n)- The 
reverse inclusion is obvious as q is a probability distribution. We now move to consider the second equality in (|10l) . 
We will here work at the level of the algebra ^4(A)® 2 , while in Appendix A the analysis is performed at the swap 
algebra C7n level. 

Let us begin by considering just two maps TZq 1 and 7ln 2 - Since one has that Fix(7£nJ = C{l,To i }(8>-4(^f)® 2 , (i = 
1,2), it is easy to check that there are two cases: 1) If fii n fi 2 = then F = C{l,T ni } (g) C{l,Th 2 } ® A; and 2) 
If fii n fi 2 ^ then P = C{l,T niU n 2 } ® A where P = Fixf^nJ n Fix(P.Q 2 ) and A := .A((fii U fi 2 ) c )® 2 . Case 1) is 
obvious, we will then consider 2). 

Let Ai, A 2 and A 3 be as in ii) in the proof of Prop. 1 (Uf =1 A; = fi x U fi 2 )). If X e FixP.^ (X e Fix7£n 2 ) one 
can write X = al x <g> 1 2 ® X 3 + /3Ti ® T 2 ® Y 3 (X = q/Xl <g> 1 2 <g> 1 3 + /3'Yi ® T 2 <g> T 3 ) where X 3 and Y" 3 (Xi and F x ) 
are arbitrary operators in *4(A 3 ) 2 ( .A(Ai)® 2 ). Therefore if X e F one must have all ® 1 2 ® X 3 + /3Ti ® T 2 ® F 3 = 
a'Xi ® 1 2 ® I3 + /?'Fl <g> P 2 <S> P 3 , this equation can be solved for arbitrary a and /3 iff li ® 1 2 ® X3 ~ Xi (g> 1 2 (g> 1 3 and 
Ti(g)T 2 (g)Y3 ~ Yi<g>T 2 <g>T 3 . This implies Xi = li and X 3 = 1 3 as well as Y x = T x and Y 3 = T 3 . In other terms the general 
elements e P must have the formX = ali(g)l 2 (g)l 3 +/3Ti(g)r 2 (g)T 3 namely X e C{l ni un 2 , 7h lU n 2 } = Fix(P ni un 2 )- 
This shows that P C Fix(7?,n lU Q 2 ), the opposite inclusion is obvious, hence P = Fix(72.f2iUQ 2 )- 

The intersection over the whole C required in (|10|) can be now perfomed, using 1) and 2), by organizing the 
hypergraph C in to its connected components Cj, (i = 1, . . . , K). Different C^s give rise, thanks to 1) to the different 
tensor factors in (flOl) ; whereas intersections PlneC Fix(Pa) give rise, thanks to 2), to Fix(7£c0 = Tc;}, where 

Cj := Uoec ^- ^he ^ as ^ factor m (HOI) simply reflects the fact that all the elements in the algebra over the set 
(Ui=i C«) c are (trivial) fixed points for all the TZq's. (Notice that often this set will be empty). This completes the 
proof of (HD). 

ii) Since 7Z is a non-negative operator on the Hilbert-Schmidt space .4(A)® 2 with unit norm its spectrum 
Sp(P) = { Pa } a is contained in [0, 1]. Therefore K k (T n ) = J2 a P k a F a {F a , T a ), (n(F a ) = Pa F a ) and lim^oo K k (T n ) = 
Eq P =1 Fa(Fa, Tq). The relevant eigenoperators F a can be chosen from an orthonormal basis of the space Fix(P.) 
that we characterized in i). This latter eigenspace, according to i) has dimension at least 2 K but symmetry further 
simplifies the computation of P^. Notice indeed that, for a totally factorized state uj = ^ieA^iji't'Ai one nas that 1) 
(w® 2 ,T ) = 1, (Vfi C A) and 2) (w® 2 ,X 1 X 2 ) = (cjf 2 , X x ) (w® 2 , X 2 ) if X x and X 2 are supported on disjoint subsets 
ill and fi 2 and u>i — l S>jecii\4'j}{4'j I) (* = 1,2). These two facts implies in turn that, in each of the two-dimensional 
C{1, Tci} factors of Fix(P.), only the symmetric term 1 + Tc i will have give a non vanishing contribution to purity 
(as (oj® 2 , 1 — Tc t ) — 0.) Equation (|TT|) then follows by computing the Hilbert-Schmidt scalar product of the swap 
To with the normalized n^i( a + T c0- From Q = (n n (uf =1 Ci) c ) u(Un (ll^Ci)) = (fin (U^iCi) ) U (u^^,) 
it follows that Tq = (^[\f = Po 4 J Psin(u K 1 c i ) c - Since 7?. has no non-trivial action on the last swap in Tq_ that term (if 
present) gives rise to an irrelevant factor one in P^. The other K swaps instead give 

p °° = n^ M ic^iSi +1) ' t *)(" 93 > V 2di4diS+i) > = n< d ic i i ( ^c i r+i)' 

computing explicitly the scalar products above (use Tr(T Ui ) = d?^-^ d) a <\ , Tr(T Ci T Qi ) = Tr(Tb i _ n . i ) = 
dl^l-l^ld 2 !^!) proves dHJ). □ 

Remark.— As noticed in the above for a totally factorized state w = ®ieA|^i)(^i| one nas that (w® 2 , Tq) = 1, VO 
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regardless of the choice of the |</>i)'s. It follows that in purity calculations u>® 2 can be replaced by the average 

iSA iSA yJ ai,...,Q| A |=0,l u + 

where d + := d(d + l)/2 and Il + := 2~l A l^ 0cA Th = |7a| _1 StgTa ^ From this l as t expression is clear that II + 
is the projector over the identity irrep of Ta acting on .4(A)® 2 . Moreover, ||n+||| = TrII+ = 2~I A I E ncA TrT n = 

2" |A| E^cA d|f2|d2(|A| " |n|) = 2 " |A ' d2|A| ( 1 + ^ 1 ) |A| ^ H n+ ll2 -rf+' /2 - 

Example .— Let Q = (A, E) be a connected graph and the ensemble map given by: 1Z = {E^ 1 J^eeE^-e- This is 
the so-called random edge model considered in [3, [l(| ■ Here we have that £ coincides with E (regarded as a family 

of subsets {vi,V2} C A) therefore one has K = 1 and \C\\ — |A| and Eq (fTTj) becomes P^ = ' ■ 



B. Area Laws for the uncorrelated case 

In this section we exploit the locality structure of the ensemble maps and the swap algebra relation ^ to prove short- 
time upper bounds for the average purity of a local region A. For simplicity, we will consider a uniform distribution q 
over £. 

From Eq. ([9]) it follows that the maps 1Z act on C7a according the rule 1Z(Ta) = J2bcA Rb.aTb, {A C A) where 

Sb,a-q + ot-(A, O) 6b,auq) ■ (12) 

necn(dAy nedA 

For a totally factorized initial state, one has Pk = (uj® 2 , 7Z k (Ta)) — YIbi B k ^-B k ,B k - 1 RB k - 1 ,B k - 2 ' ' ' Rb 1: a- If 
C := |Q n A\ = |fi - A| then c~ = =► a+ = a_ = c+/2 = d c /(d 2C + 1). In the rest of the section we will 
assume that a- = a+ =: A^, (VT2 E £, A C A) [25(. The matrix (p~2|) can be written as a sum of a diagonal and a 
off-diagonal matrix R = F + Fi where F := diag (1 -p{A)) A cA, (p(A) := EnedA ^(P) = 1 - Y,necn(dAy and 
(^i)b.a = Nd J^nedA ?(^) ($B,A-fi + <5b.Auo) • The off-diagonal matrix i*\ connects, with strength iVj, each region 
A with other 2|cL4|. For example, for k — 1 one immediately finds 

p i = ^ = H (( F o)b,a + {Fi) b ,a) = 1 - P{A) + 2N d = 1 - P(^)(l - 2 ^z), 

bca bca neo^t 



B.A 



For a uniform distribution g over £ one has p(A) = |cM|/|£|. In general: Pk = J2bcA a k =o ^ai ' ' ' ^q* 

E[a]ez| Usin S the matrix structure ([11]), one obtains Y { A k) {a) = ^ BcA (f Ql • ■ ■ F ak ) B ,A < (2N d p(X))\ a \(l - 

p(X)) k -^ where: \a\ := £*L X a h x ( X ) is the region in the family obtained by A with k successive allowed operations 
[26j given by either joining or subtracting elements of £, with the largest (smallest) \dX\ (\dX\)[27\ Whence 

fc k 
ft < V C fcjQ |(2iV d p(X))l«l(l-p(A > )) fc - |a| = (l-p{X)+2N dP {X)) =(l-(l-2N d )p(X)+p(X)-p(X)) k 



\a\=0 



= (1 - e pP (X)) k (l + AkA (Y , ) k < exp 
1 - e n p(X ) 



k(p(X)e p -p^) 



e p p(X) 

where i a i are binomial coefficients, = p(X) — p(X) and 



(13) 



I'd - l) 2 

e p ;=l-2N d = { -^ T L. (14) 

This constant is the Haar average of the "entangling power" for unitaries acting on d x d— dimensional systems 
introduced in [19(. If we assume that each step the boundary size of a region can change just by a quantity 0(1) = 
o(\dA\) i.e., the hyper-graph £ has bounded degree and the boundary size of A is large scaling quantity, it follows 
A fc = 0(fe/|£|) whereas p{X) = 0(\dA\/\C\). In words: inequality (Til?l) describes, for k = O(l), an exponentially 
bounded decay of average purity with a rate that fufills an area law (up to corrections 0(1)) and that is proportional 
to the average entangling power e p 
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C. An exactly solvable model: the uniform path-graph 

In this section we wil discuss a one-dimensional uncorrelated case where the use of swap algebra formalism allows for 
an explicit exponential reduction of computational complexity [28| . The ensemble map TZ can be explicitly diagonalized 
in the relevant invariant sub-space of C7a and bounds on the convergence rate Pk — > Poo ar e easily established. 

Let G = (A,E) be a path-graph of length L : A := {1, . . . , L}, E = {{1,2}, {2, 3},...,{L - 1, L}}. From Eq © 
it follows that the (L + l)-dimensional subspace of C7a spanned by {T^ / 51 = {1, ...,£}, i = 1, ... L} and T$ = 1 is 
invariant under all the edge maps TZ e (e £ E). We denote the basis of this (L + l)-dimensional invariant subspace by 
{|0), |1), |2), . . . , \L)\}. In the uniform case one has TZ = J2i=i an d relations (fT2f give: 

TZ\0) = |0>, TZ\L) = \L), Tl\i) = o|i) + 6(|t-l) + |t + l)), (i = - 1) (15) 

where a := (L — 2)/(L — 1) = 1 — (L — 1) _1 and 6 = AT^L — 1) _1 . Let us denote by R the (non-hermitean) 
matrix representation of TZ in the basis {|i)}^_ . Eq (|15| implies that P is a sum of two matrices i.e., R = R± + R% 
where Ri = |0> <0| + fe|0)(l| + \L)(L\ + b\L)(L - 1| and R 2 = al L _i + bA t -i. Here l L _i denotes the identity 
matrix in the (L — 1)— dimensional space spanned by {li)}^ 1 and Al_i is the adjacency matrix of the path graph 
of length L — 1 whose nodes are labelled by the same basis vectors. The matrix R isdiagonalizable by a (non- 
unitary) transformation R = P~ 1 DP where D = diag(A^)^_ an d P is a, matrix whose columns are normalized 
eigenvectors of R, i.e., P jtfl = = A p |*) M , (j, = 0,...,L) [29] The spectrum of P is: {A M } M := 

{1, a + 2&cos( 2 j i ), (/i = 1, . . . ,L — 1), 1}. The L — 1 eigenvalues different from one can be written as 1 — where 
A h := 2^(1 - 27V d cos(2^)) e (0, 1). The spectral gap A := 1 - max{A ft }^ 1 1 is then given by 

A = minAfc = (l - 2JV d cos(^)) > ^(1 - 2N d ) = ^. (16) 

Here e p is the average entangling power defined in Eq. . The lower bound (|16p shows that the spectral gap of the 
ensemble map TZ of the unifom path graph is "large" in the sense that it does not vanish faster than for large 
system size L. This fact, in turn implies a fast convergence rate Pk i— > Poo for k — > oo. 

For general k and initial state \l) one has P fe = Y,j=o( Rk )j,l = T,j=o( PDP ~ 1 )j,l = T^=o{ PDkp ~ l )j,l = 
E^o^E^o^c)^ Since = a l = 1 and |AJ < 1 for /z = 1, ...,£. — 1, one obtains P k - Poo = 
Y^hZx (1 — ^ft) fe (Sj'=o p j.h){ p ~ 1 )h,i- A direct computation shows that 

2 v — > , ,i. ,nlh . 

P k - P oo = -j- 22 t 1 - A '') Sin ( — ) 
h odd 

where the sum is over the odd numbers h between 1 and L—l, and, consistently with (fTTj) , Poo = (d L ~ l +d l )(d L + 1) _1 
[2(ij |. From this equation and the spectral gap definition (ITBI a straightwforward argument shows that |P/- — Poo| < 
le~ kA C where C = O(l) is a constant depending on <i. Then, using the spectral gap lower bound in (TTBl . one obtains 
the convergence rate estimate 

k> — (log(-)+ log 0=> |P fc - Pool <e. (18) 
e p e 

This inequality shows that, for fixed accuracy e, the circuit length k scales linearly with the total system size [Tl| . 
Interestingly (|18[) also shows that the convergence rate is pro portional to the the average entangling power e p , an 
intuitive result given the physical meaning of this quantity [19(. 

We would like to end this section by showing that in this one-dimensional case the bound (fT3"| (with Afe^ = 0) is tight 
for for short times k. To begin with, notice that if k < min{Z, L — l} the action of the map TZ k on \l) does not involve the 
boundary states |0) and \L). As long as these latter fixed states can be discarded Eq. (flip)) shows that TZ acts as the sum 
of three commuting translation operators i.e., TZ = 5E a =o ±i c aT a where T a \l) = \l + a), (I = 1, . . . , L — 1; a = 0, ±1) 

and Co = a, c±\ — b. Therefore TZ k \l) = J2 ai a k =o ±i c «i ' ' ' c a k T ai ^ \-a k \l)- As remarked in the above, for a totally 

factorized initial state, each vector in the last formula gives a unit contribution to the average purity Pk, hence 
p k = J2 ai a k =o ±i Cq i ' ' ' Ca k = ( c o + c -i + c i) k — ( a + 26) fe = (1 — j^e p ) k . This simple result can be also checked 
by resorting to the explicit formula (IT71) . 



2NdM?£) nh, 

-A h cotf — , 

2jV rfC os(^) - 1 K 2L' 



(17) 
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D. Correlated case 



In this subsection we will consider an ensemble map of the form 7Z = Yinec This is a natural generalization of 
the correlated family of 1-dimensional L-RQC called the "contiguous edge model" introduced and analyzed in [9l [i"o|. 
In the contiguous edge model the graph vertex set is given by A = { — L, —L + 1, . . . , —1, 0,1..., L} and the edge 
set E as in Sect. C. The L-RQC model is now defined as follows: At each discrete time t = k one selects 2L (Haar) 
random U e and U is built according to a permutation ir e $\e\i = 2L) U = U e7r{1) U e7r(2) ■ ■ ■ C/ e7r(|E|) and then 
the process is iterated at t = k + 1 and so on. Each U e . is unitary acting on the edge ej — (J, j + 1) state space 
hj ® hj+i = (C d )® 2 . The ensemble map corresponding to the specified model is given by TZ n :— 1Z e ^.. En o • • • 72-e„ (1) ■ 
Concerning the map fixed point, in [lj| it was conjectured that for all n e S\e\ the relevant fixed point of TZ^ is given 
by 1 + Ta and Poo by (|11[) . In the following we will demonstrate a result (Prop. 5) that, as a particular case, proves 
these conjectures. 

Notice that, since in general different TZq do not commute, the map 1Z = Ylnec ^ is not even necessarily 
Hermitean. Therefore one needs a different type of analysis respect to the uncorrelated case of Prop. 4. 



Lemma 

Let Pi, . . . , P„ be projectors in a finite-dimensional Hilbert Space, R —: P„P„_i ■ ■ ■ P\ and F := n" =1 ImPi. Then 
i) Fix(P) = P, ii) ||P|| = 1 <S> Fix(P) ^ {0}, iii) R = P + Q where P is the projection on P and ||Q|| < 1. 

Proof.— i) Obviously P C Fix(P), let us show the reverse inclusion. Let xo € Fix(P) and Xi — P^Xj-i, (i = 
l,...,n). One has x = R(x ) = x n = P n (a; n _i) and then ||a;o|| = ||P„a;„_i|| < ||P„|| ||a;„-i|| < ||P n ||||P T ,_i||||x 7l _2|| < 
"* — n"=i ll^llll^oll < ll^oll: this shows that the equalities ||PjXj_i|| = ||xi_i||, (i — 1, . . . ,n) hold. These latter, since 
the P^s are projections, imply Xi :— PiXi-i = Xi-\ (i = 1, . . . , n). Therefore x n = x n _i = ... = x\ = x$, namely 
PiXQ = xq for all i i.e., xq g P. 

ii) If P ^ {0} obviously ||P|| = 1; let us show the reverse implication. If ||P|| = 1 from compactness of the 
unit sphere 3x* , ||x*|| = 1 s.t. ||P.t*|| = ||a;*||. From this relation, following the same steps as in i), one sees that 
P(x*) = x* for all i = 1, . . . , n i.e., 0/i'ef. Since ||P|| < n"= ll^ll = 1, if P = {0} then ||P|| < 1. 

iii) P is a subspace of all the spaces ImP^, (i = 1, . . . , n) hence one can write p = P + Qi where Qi projects on the 
orthogonal complement (in ImP; )of P. Now we have R = (P + Q n )(P + Qn-i) ■ ■ • (P + Q\) = P + Q n Qn-i ■ • • Qi ='■ 
P + Q where we have used QiP = PQi = 0, (i = 1, . . . , n). Notice now that Q = R(l — P) therefore Q is a product of 
projections that (by construction) cannot have a non-trivial joint fixed eigenspace. From ii) it follows that 3A £ (0, 1) 
s.t. IIQH =1- A. " □ 



Proposition 5 

i) Let 1Z — 72.n|£|7^fi|£|_i • ■ - TZni, ({"-i/iJi := £) be the ensemble map of a correlated L-RQC. For a pure totally 
factorized initial state ui one has 

P fc = (^ 2 ,^ fe (Po))-Poo+0(2 |A|/2 e - Afc ), A:-lo gT ^, (19) 

where is given by pT|) and 1 — A = \\1Z(1 — V)\\ < 1, V projector on Fix(P.) (see Prop. 4). 

ii) If IZ-k := "^•n„n£i)'^'^ W (|£i-i) ' ' ' ^ n 7r(i) nas the same infinite time purity Eq. (fl9)) \/ir £ S\c\- 

Proof.— i) We are in the setting of the Lemma with P 7^ {0} (Prop. 4) where now the finite-dimensional Hilbert 
space is tye Hilbcrt-Schmidt space £(A)® 2 the projections are the TZq and n — \C\. 

One hasPJ 5 = {V + Q) k = V + Q k {TQ = QT = 0). Since ||Q fc || < ||Q|| fc = (1 - A) k = cxp (-klogj^j it follows 

that limfc^oo Q k = and hence lim fc _ ! . 00 lZ k = V. Therefore P k = (w® 2 , (V + Q k )(T n )) = P x + (cj® 2 , Q k (T n )), where 
Poo = (u;® 2 , P(Psi)) is given (see Prop. 0) by Eq. (|11[) . Now, using the Remark after Prop. 4, and Schwarz inequality 
for the Hilbert-Schmidt scalar product we can write (w® 2 ,Q fc (T n )) = d+ |A| (n+, Q k (T n )) < d~ lM \\U + \\ 2 \\ Q k (T n )\\ 2 < 
d -|A|/2|| gfc(Tfi) || 2 < d;l A l/ 2 || Q fe||||T || 2 < (^)|A|/2|| Q ||fc < (-2^)IM/2(i _ A ) fe < 2l A l/ 2 (l - A) fe ). Where we also 
used HThll 2 , = TrP^ = Tr 1 A = d 2 ^ . This proves (fT9|) . 

ii) Point i) shows that infinite time purity Poo depends just on the joint fixed space n'£' 1 Im7^ 7r (i) = ImP. This 
latter is clearly independent on the permutation n. □ 



In the decomposition TZ^ =V + Qk the second term does in general depend on ir Accordingly the convergence rate 
of P. to Px, (related, by Prop. 6, to HQ^H) may depend on tt. This phenomenon has been numerically observed in 
@,Gl}. From ((Til) it follows that k > A -1 log(2' A ' / ' 2 /e) => \P k - P^] = 0(e). This shows how the convergenge rate 
Pfc — > Poo is controlled by the norm ||Q|| = 1 — A. This norm is also the second largest singular value of 1Z. The 
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positive constant A is therefore just the difference between the two largest singular values of the ensemble map 7Z 
and it is often referred to as the spectral gap (see Sect C). From the practical point of view it is crucial to determine 
how the convergence varies with the ype of L-RQC and how the spectral gap scales as function of |A| [ill, EH- I n 
and [Io| are given numerical results for a few natural choices of n. In particular for the so called expanding n (see 
Sect V.A in |l0(), arguably the less efficient at generating entanglement, it was found P k = 2(1 — e p ) fe (l + e p )~ h (e p 
is defined in ([M)) for 1 <C fc < |f2| and a cross over to a volume- law for k = Moreover: 1) asymptotic purity is 
given by (TlT|) (for K — 1) and, 2) the spectral gap of TZ^, was found (numerically) to have a finite limit, 1 — (2Nd) 2 , 
for L — ¥ oo. 

We end this section by noticing that, for both the uncorrelated and correlated cases, Propositions 4 and 5 show 
that the fixed point problem for the ensemble maps 1Z can be formulated as a ground state problem for the ensemble 
"Hamiltonians" K, := Ylnec^ ~ ^fi)* This are local non-negative operators acting on the Hilbert-Schmidt space 
A(A)® 2 = <g) igA 2?(/if 2 ). In fact one has /C$ = 71q& = $,Vfi e £, namely the ground state manifold of JC 
coincides with the 2 x -dimensional space Fix(lZ). It would then interesting to see whether one may adapt to the 
present context the techniques developed in [2l| to lower bound spectral gap of local Hamiltonians as well as to see 
whether one may exploit in convergence rate estimates the so-called quantum detectability lemma (22j . 



VI. ASYMPTOTIC STATISTICS AND LOCAL STATE DESIGNS 



In this section we will establish a sim ple connection between the formalism and the results of this paper and the 
theory of quantum local t— designs 0, HUH- 

The asymptotic average purity value is 0(e?~'° ) away from the minimal possible one for the reduced state of 
the region f2 i.e., P m in = gH 51 '- This implies that for sufficiently large k most of the states wn(Ufc) := tr^c(UfcU;Ut) 
(Ufe RQC of length k) should be close to the maximally mixed state ln/d' n '. Indeed on can prove that if k is such 
that \P k - Pool < e then 

IM U *) - ^| h" < V^/^+- (20) 

Therefore if \Q C \ = (1 + a)\Q\ (a > 0) and e = gH° ' one has that the average trace norm distance between states 
generated by the ensemble of RQCs and the maximally mixed one is 0(d" a ^ n ^ 2 ). Using Markov inequality, for large 
|fi|, one sees that with high probability each state in the ensemble is locally nearly indistinguishable from the totally 

mixed one. Notice now that similar claims can be made for the Haar measure as: trwji(V-') 2 ^ < ^[qj [HI where 

All this suggests that statistics of local observables in fi over the L-RQC ensemble (for large k) should be a good 

approximation of the Haar's one. In order to see this one has to show that the distance T> t := Hw^Ufc)®' * — 

wn(V , ) (8t ||i; (i-positive integer) is small e.g., less than S in some limit. If this is the case all local observables will 
have all the L-RQC statistical moments up to the t-th one at most S away from the corresponding Haar's one [?}■ 
An ensemble of states having this property i.e., to provide a way to approximately sample the Haar distribution is 
called an (^-approximate state t-design. How this can be done by efficiently resorting just to local quantum circuits is 
a problem currently at the centre of an intense attention in view of its applications to quantum information theory 
(see recent [§]). 

Proposition 6 For k = O (A -1 logd(|A| + (1 + a)|0|)) , (a > 0) the L-RQC family is a S— approximate state 

t-design for local observables A € A(Q) with 5 — 0(V te~ a \ n \). 

Proof.- For k = O (A" 1 logd(|A| + (1 + a)\Q\)) bound (??) implies \P k - P^l < := e Using triangular 

inequality, convexity of the trace norm, Eq. ([20)) and — < Vt\\<-o — <j\\i (for all pairs of density ma- 

trices p and a) one obtains V t < ||w n (U fc )® *- Jft® ' ||i + 11^®* - w fi (^)®* ||i < ||w n (U fc )®* - Jft 8 *||i + 
\\^ m -^n < V~t + dMe + yf^j =: 5 = O(v^W) □. 
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VII. CONCLUSIONS 



In this paper we provided a mathematical presentation of the approach to local random circuits discussed in 
, [13] ■ We introduced different classes of local random quantum circuits and showed that their statistical properties 
are encoded into an associated family of completely positive maps. Ensemble average of quantum expectation can 
be studied as a function of the circuit length and we proved infinite time results as well as convergence bounds. 
Remarkably, average entanglement dynamics can be described by the action of ensemble maps on operator algebras 
of permutations (swap algebras). Also in this case we proved infinite time results for the expectation value of the 
purity of a local region both for uncorrelated and correlated local random quantum circuits and short time area-law 
bounds for the uncorrelated case. The swap algebra formalism is poweful as it allows in some case an exponential 
reduction of the problem size i.e., d 2L k> L. We illustrated such a phenomenon by an exactly solvable one dimensional 
uncorrelated case. Finally we briefly discussed a connection with t— design @, H[ for localized observables. 
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[23] For example: fli C fl 2 => Imft^ D Imftn 2 TZn 2 (1 - 7£nJ = 0. 

[24] A connected component in £ is a maximal set d C Ung^fi C A such that Wu,v £ Ci, 3fii, . . . , fin £ £ such that 
u £ fli, u 6 fin and fij n Qi+i ^ 0, (i = 1, . . . , n — 1) i.e., given any two vertices u and v there is a path of hyperlinks 
connecting them. 

[25] This is the case when the elements fii in C are edges of a graph Qi — {vi , v 2 } (vi, v% £ A). In fact fix £ dfl => f2i — f2 = {vi}, 

and Slifl!! = {« 2 }. Here N d = d/(d 2 + 1)0. 
[26] By allowed we mean that at each stage the element fl £ £ added or subtracted to A' has to belong to dA' 
[27] Notice that X and A may depend on both A and k. 

[28] In 10] it was analyzed the case of the complete graph with L vertices. Also in this case one has an exponential reduction 
of complexity as average purity dynamics can be restricted to the (L + 1)- dimensional S^-symmetric subspace of C7a — 
(C 2 )® L . 

[29] Explicty one finds: |¥ ) = |0>, \V L ) = \L) and |* h ) = N h (^=1 Ch0')li> ~ b /A h (Ch(l)|0) + c h (L - 1)|L))) , 
where Ch(j) := sin(nhj/L) for /i = 1, ...,L — 1. The normalization constant A/h is given by A4 = 
(L/2 + (b/A h ) 2 (cl(i) + c h (L - 1) 2 )- 1/2 = 0(L" 1//2 ). 



Appendix A: Fixed swaps subalgebras 

In this section we will show how some of the facts used in the Proof, of Proposition. 4 can be established working 
entirely within the swap algebra formalism. 

Proposition Al 

Le TZq 1 be an ensemble map given by ©, Ox , 2 C A and Fix (TZq^ := span{Tn / TZq 1 (Tq) = Tq}. 

a) Ai C A implies that C7ai = span{Tb / C Ai C A} is a subalgebra of C7a of dimension 2' Al L 

b) if Ax n a 2 = cr AlU A 2 = cr Al v cta 2 = cr Al ® cr A2 . 

c) Fix(^J = span {T n /OD Oi V O C Of } = span{T n /Oi ^ 90}, dim Fix {TZqJ = 2l A H^il+i 

d) Fix(7£m) is a subalgebra of C7a isomorphic to C{1, Tq x } €5 C7nj. 

e) n x n 2 = Fix(7£ Ql ) n Fix(^o 2 ) S C{1, T n J ® C{1, T n J ® CT^un,)* 

f) Hi n 2 ^ Fix(7e ni ) n Fix(^n 2 ) = Fa^un,) S C{1, T niU n 2 } ® C^un,)- 

Proof. 

a) Obvious. 

b) Any element of Ai U A2, when Ai n A 2 = 0, can be written in unique way as Oi U 2 where 0j C Aj (i = 1, 2). 
Therefore Tn lU n 2 = T ni An 2 = Tq x Tq 2 showing the first isomorphism. Moreover, the map l{Tq iU q 2 ) = Tq % ® Tq 2 
can be checked to provide the second one. [Notice also dimCTAiUA;, = 2' AlUA2 l = 2' Al 2' A2 ' = dimCTAi dimC7A 2 -] 
If Qi,fii £ Ai(i = 1,2) one has (Oi U 2 )A(0i U 2 ) = (O1AO1) U (0 2 A0 2 ); it follows that t(T 0lU n 2 ^unJ = 
t (^'(niuo 2 )A(o 1 u^ 2 )) = t ( T '(ii 1 Ao 1 )u(n 2 Asi 2 )) = ^niAfii ® r o 2 Ao 2 = Cfoi ^ ^ 2 ) (T^ i ®T^ 2 ) = t(To lU f2 2 ) t(T^ iLl fj 2 ) i.e., 
t is also an algebra isomorphism (not just a vector space one). 

c) Notice first that, in view of © , Fix (TZq 1 ) is a linear subspace of CTa that contains span {Tq /fiDfliVOc Of} 
(as for every element one has Oi ^ <90). Moreover this latter subspace has dimension 2l A H 0l l +1 as every element 
of its basis has the form 0' or 0' U Oi where 0' C Of. Since 1Zq 1 is a projector the dimension of Fix (TZq^ i.e., the 
eigenspace of TZq 1 with eigenvalue one, is given by TtTZq 1 i using again © one sees that this trace is in fact 2l A H 0l l +1 
(only the elements such that Hi ^ dft contribute, each by a one). This proves the equalities in c). 

d) If fl, fi' are such that 7b, 7b< £ Fix (TZq 1 ) there are four possibilities 1) fl C 0? A fl' C fif 2) fl D Or A fl' C fif , 
3) il C f2£ A f2' D f2i, 4) i7 D A fi' D f2i. In all these cases it can be directly checked that flAH,' either 
contains fli or is contained in il^ i.e., f2i ^ 9(51AJ7'). Hence TqTq' = Toao' G Fix(72.ni)- The isomorphism in d) 
follows from the fact (noticed in c)) that every element in Fix(7?-o 1 ) can be written as Tn^w = i (Tq 1 (g) Tqi) or 
T iV = L A (tSiT sr ), (Q'cn'i)- 

e) If O e Fix(ft 0l )nFix(fto 2 ) (with O x n 0. 2 = 0) there are four cases: 1) D Oi AO D fi 2 , 2) C AO D 2 , 
3) HdHiASIcHjJ) ft C Of A O C 2 . Is easy to see that these alternatives account for the four different types of 
terms one gets from the tensor product in e) e.g, i~ 1 (Tn 1 ® Tq 2 ig) Tq>) — Tq iU q 2U q>, (O' C (Oi U 2 ) c ) corresponds 
to 1) and so on. 

f ) If OiD0 2 ^ then cases 2) and 3) in the above are not allowed and one is left with 1) and 4); namely, O D OiU0 2 
or O C Of n 2 = (Oi U 2 ) c . This means Oi U 2 ^ 90 that, thanks to c), is equivalent to say Tq £ Fix(7?.si lU o 2 ). 
Last isomorphism follows from d). □ 



